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ABSTRACT 

The nuclear and emitted radiation characteristic a of the radioisotope 
elements and impurities in commercial grade plutonium dioxide are 
presented, in detail. The development of the methods of analysis are 
presented. 

Radioisotope thermoelectric generators (RTG) of 1575, 34.68 and 5679 
thermal watts are characterized with respect to neutron and gamma photon 
source strength as well as spatial and number flux distribution. The 
results are presented, as a function of detector position and light element 
contamination concentration for fuel age ranging from 'fresh' to 18 years. 
The data may be used to obtain results for given and ^®Pu con- 
centrations . 

The neutron and gamma photon flux and dose calculations compare favorably 
with reported experimental values for SNAP-27 
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1 . INTRODUCTION 


The ever increasing requirement for long-lived and compact power sources 
has led to the recent evolution of a number of very successful generators. 
Specifically , the development of the radioisotope thermoelectric generator 
(RTG) for both terrestrial and space flight applications has progressed very 
significantly. A number of plutonium dioxide fuelled RTG's are already in 
actual use. The NIMBUS weather satellite is powered by a Space Nuclear 
Auxilary Power system, i.e. , SNAP 19. Instrumentation left on the Moon 
by the Apollo 12 and 14 missions are powered by SNAP 27 units. Future 
space missions include RTG assemblies. Unlike solar cell generators 
their power output is not dependent upon solar flux and thus they are 
suited for deep outer planet space missions of 10 to 20 year duration. 

The radioactive properties of RTG's require that special attention be given 
to their engineering design from both the standpoint of mission safety and 
mission operational requirements. This report is concerned with these 
factors as they relate to the gamma photon and neutron radiation emitted 
by plutonium -oxide fuelled generators. The radiation data presented will 
be of .interest to spacecraft, RTG and science experiment designers. It 
permits radiation field maps to be developed to aid in location and protec- 
tion of sensitive electronic packages and science experiments . The data 
may also be used for determination of personnel exposure and protection. 

The RTG gamma photon and neutron radiation result from the natural radio- 
activity of the plutonium dioxide heat source. The gamma photons are due 
mainly to the plutonium isotopes and their daughter products with small 
contributions resulting from induced fission and neutron interactions with 
materials in the source assembly. The neutrons result mainly from 
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plutonium decay alpha emission. and thus (a, n) reactions with low atomic 
number inpurity elements including the O component of the oxide; a 
small contribution is due to plutonium spontaneous fission. The neutron 
flux may be decreased by reduction of the impurities and depletion of 


The natural radioactivity of the isotope 238 Pu consists of alpha particles 

accompanied by a spectrum of gamma photons and spontaneous fission 

neutrons. Commercial grade plutonium consists of 238 Pu and other Pu 

isotopes as indicated in Table 1-I^" 3) . Even though the 238 Pu is only 

~ 80% abundant, it contributes ~ 99.9% of the total alpha activity. Since 

the average energy of emitted alphas is ~ 5.5 MeV, (a, n) reactions occur 

with low atomic number elements such as oxygen which is abundantly 

present in Pu0 2 as indicated in Table 1-1, The use of oxygen in which the * 8 0 

isotope is depleted by a factor of 100 reduces this effect correspondingly by~ 100. 

Again, the presence of impurity beryllium or fluorine even in low anundance may give 

rise to a significant (a,n) neutron yield; some light impurity element yields are 

(3) 

given in Table l-II . The spontaneous fission (SF) neutron emission 

238 

of plutonium is primarily that of Pu since the other contributing isotope 

240 , 23fl 

Pu, has an SF rate ~ 0,015 that of Pu. The SF neutrons have a 

Maxwellian distribution peaked at~0.7 MeV, an average energy of ~ 2.0 MeV, 

and a maximum energy of ~10 MeV. Fast photo -neutrons are also pro- 

2 3ft 

duced in PuO . due to energetic gamma photons from Pu spontaneous 
“226 208 

fission and the Pu daughter, T^, interacting with impurity 
elements . 

Pu0 2 sources emit X-rays and gamma photons as a result of spontaneous 
and induced fission, non-fission neutron interactions, decay of fission 
products and their daughter products , and alpha Interactions with light 
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1 8 

elements such as O. The gamma photons emitted by PuO are presented 
in detail in Section 2 of this report. The gamma photon activity from Pb, 

O QO 

Bi and Tt daughters increases with time after Pu purification to a 

maximum at ~ 18 years; this is readily seen in the decay in 

Figure 1-1, Freshly separated plutonium fuel has a 'soft' gamma 

spectrum. With increasing time the decay daughter products of the small 
236 ^ 

Pu impurity increase the intensity of the 0.2 - 0.3 and 2.0 - 3.0 MeV 
photon groups. Spectral hardening is complete at approximately 18 years. 

241 241 

Long-lived Am, a daughter of 13 year Pu, yields increasingly intense 
gamma photon emission rates with age. However, since the photon energy 
is 0,06 MeV, it is easily shielded and so is not intensely present in the 
emission spectra of encapsulated sources. 

As indicated above fast neutrons are emitted by PuC> 2 primarily as the 
result of the three interaction phenomena: spontaneous fission, alpha- 

fa) 

neutron emission rates for these effects; detailed neutron data is pre- 
sented in Section 3 of this report. The (a, n) neutron emission rate depends 
on the impurities present in the source, since the reaction requires the 
presence of low Z elements, as seen from Table 1— II. Table 1 —III lists 
neutron activity due to (a, n) for normal and O depleted PuO.; it also 
gives the photo-neutron production rate. 

The alpha and beta radiation emitted by PuOg sources does not 
present any problems if the source encapsulation has a wall thickness 
greater than ~ 1000 mg/cm^. (e.g., ~ 0.06" Ta). 

Section 2 of this report presents gamma photon source and calculated 
emission flux data for a hypothesized SNAP RTG. The data was deter- 
mined for 1575, 3468 and 5679 watts of thermal power. The flux data 
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was calculated as a function of PuO^ age at the RTG exterior surface 
and points on both the axial and radial mid-planes, Section 3 of this 
report presents neutron source and emission flux data similar to the 
photon data of Section 2. Section 4 consists of a brief review of the data 
in Sections 2 and 3. A comparison of the calculated data with published 
data is presented. Pertinent information and methods of calculation are 
given in Appendices A through J 
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TABLE I-I 


RADIATION PROPERTIES OP PuO 


(1-3) 


A. 


Chemical Composition PuO 


2 


Plutonium (total) 88,2% 

Oxygen 11 , 8 % 


B. 

Isotope 

Isotopic Composition^ 

Activation of 
, Pure Isotope,. 

Abundance %= (a) d/sec-g = (b) 

Contribution to the 
Activity of 1 gram of 
Product® (a) x (b)/l 00 

236 Pu 

1.2 x 10" 4 

1,97 x 10 13 

2.37 x 10 7 

(a) 

238 Pu 

81 % 

6.35 x 10 U 

5.14 x 10 11 

(a) 

239 Pu 

15 % 

2,27 x 10 9 

3.41 x 10 8 

(a) 

240 Pu 

2.9 % 

8.38 x 10 9 

2.43 x 10 8 

(a) 

241 Pu 

0.8 % 

4.12 x 10 12 

3.30 x 10 10 

(0) 

242 Pu 

0.1 % 

1.44 x 10 8 

1 .44 x 10 8 

(«) 

C. Specific Activity 

238 Pu 16.8 Curies/gm 


D. Density of PuOg 

3 

10.0 gm/ cm 




Plutonium Fraction 8.82 gm/cm' 


( 238 Pu 7,15 gm/cm 3 PUO 2 ) 


Oxygen 


1. 18 gm/cm" 


E. Activity of Fuel 

7.15 gm/cm 3 x 16.8 Ci/gm = 


2 3 

1.2 x 10 Curies/cm -PuO, 
4. 5 x 10 * 2 a/cm 3 - sec 


= 3.92 watts/cm" 


TABLE I- II 


SPECIFIC NEUTRON YIELDS FROM 
LIGHT ELEMENT IMPURITIES^ 


Element 

Li 

Be 

B 

C 

N 

O 

F 

Na 

Mg 

A1 

Si 

P 

s 


Neutrons Per Second* 
for one Part per Million 

4.6 

133 

41 

0.2 

0.0 (o', n threshold too high) 
0.1 
18 
2.2 
2.1 
1.0 
0.2 

< 0.03 

< 0.03 


* 


In 1 gm Pu metal containing 81% 


238 


Pu 
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TABLE Mil 

PuC> 2 SOURCE NEUTRON ACTIVITY* 3 * 

(i) Spontaneous Fission Rate 

238 Pu : 2.6 x IQ 3 n/see-gm 238 Pu 

fuel • 1.8 x IQ 3 n/sec-gm PuOg 

(ii) (<y, n) Neutron Emission Rate 

PuOg (normal O) * 1*13 x 10^ n/sec-gm 238 Pu 

18 238 

PuOg ( O reduced to 1/100) : '*'1..13xl0 n/sec-gm Pu 

(excludes (a, n) due to other Impurities) 

A 

(iii) Photo-Neutron Production Rate: ~5.0 x 10 n/sec-gm 

208 no c 

(varies with concentration of T-t daughter of Pu^ ) 


7 




2 3 8 _ 
Pu 


-fr 


2 . 


GAMMA PHOTON RADIATION CHARACTERISTICS 


2.1 Introduction 

This report section presents the calculation of the gamma photon flux 
emitted by plutonium dioxide fuelled RTG capsules, Gamma photon fluxes 
were determined as a function of energy and fuel' age for RTG model cap- 
sules with thermal power capacities of 1575, 3468 and 5679 w(th). The 
calculations included the determination of source data to account for 
plutonium decay, nuclear fission and the (a, n) reaction. 

Production grade plutonium dioxide fuel consists of a mixture of plutonium 
isotopes having the approximate composition shown in Table 1-1^; other 
possible compositions may be found in the references (2, 3). Each of the 
isotopes listed in Table 1-1 have the properties of natural radioactivity 
and neutron fission as given in Table 2-1^' ^-11)^ determination Q f the 
total gamma photon source intensity and spectral distribution is described 
in Section 2.2. 

The data presented in this report were determined for a standard RTG, 

namely: a hypothesized SNAP 27. The radioisotope fuel used in the SNAP 27 

fuel capsule is plutonium in the form of plutonium dioxide microspheres. The 

fuel is contained in a cylindrical annulus as shown in Figure 2-1. The basic 

model capsule was assumed to produce 1575 watts of thermal power. The 

1575 w(th) capsule power was extended volumetrically to 3468 w(th) anu 

5679 w(th) in order to study the effect of fuel self shielding on radiated 

photon fluxes. Gamma photon transport in the RTG was determined by the 

( 12 ) 

method of discrete ordinates . RTG surface fluxes were extended to 
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radial and axial midplane spatial detector locations by means of tho point 
( 1 3 ) 

kernel technique . The calculation of the photon flux distributions is 
described in Section 2.3. 

2.2 Gamma Photon Source Spectra 

* 

The gamma photons emitted by plutonium dioxide were separated into four 
components, categorized according to their origin, namely* 

23 S 

I. Plutonium isotopes and their decay products, excluding Pu. 

II. Fission and decay of fission products. 

2 36 

III. Radioactive decay of Pu and daughter nuclides. 

18 21 

IV. The alpha particle reaction, O (a, n) Ne, which yields a 

22 

gamma spectrum due to the decay of the Ne compound nucleus 

* 21»t 

to Ne. 

Other categories whose contribution to the RTG external radiation field 
were reasonably neglected are , namely: 

V. Bremsstrahlung and inelastic gamma photons which are a function 

of spacecraft geometry and materials. 

17 20 

VI. Gamma photons resulting from the O (&, n) Ne reaction. 

The gamma photon calculations were carried out for twenty (20) energy 
groups. The upper and lower limits of the highest and lowest energy 
intervals were taken as 7.0 MeV and 2.0 keV, respectively. In each of 
the corresponding gamma energy intervals , the sources of prominent gammas 
emitted by the plutonium dioxide fuel were determined for time periods 
beginning with the "fresh" (time of chemical extraction) fuel and extending 
up to 18 years of fuel aging. 
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The gamma photon spectrum emitted by 238 Pu, 239 Pu, and 24 l Pu was 

calculated as the abundance multiplied by the activity of the isotope at the 

equilibrium time considered. The activity of daughter nuclides 241 Am, 

2 37 2 3 7 * 

U and Np is calculated in Appendix A. The gamma activity of the 

236 

plutonium isotopes , exclusive of Pu, and their decay products varies 
only slightly up to 18 years and thus was assumed to remain constant 
with time, The gamma photons emitted by the plutonium isotopes and 
their decay products are given in Table 2-U (1 ' 3 ' 9 ' l4 K The reference (9) 
data was used where available, otherwise references (1) and (14) wore 
used. 

The gamma rays omitted by prompt fission wore assumed to have the spectral 
shape provided by the spontaneous fission of 235 U (15) . In the gamma photon 

energy range, 1 .30 to 1.0 MeV, the photon spectral distribution was obtained 
from the relationship 

N(E) * 26. 8e 2 * 38 ^, photons/fission. 

In the range, 1.0 to 7.0 MeV, the distribution was obtained from the • 
relationships 

N(E) = 8,03 ***^, photons/fission. 

Equilibrium fission product gamma photon emission was determined using 
analytic functions fitted to the data of reference (14). The functions and 
the reference data are given in Table 2 -III. The gamma photons emitted by 
prompt fission and equilibrium fission products were integrated numerically 
in each energy interval to yield the total number of photons produced per 
fission. Details of these calculations are given in Appendix B. 


236 

In the present work the Pu isotope was considered as a PuQ^ fuel 
impurity which ranged from l,2ppm (1,2 x lo" 6 gm 236 Pu per gm PuOg) , 
as typical of the current commercial product grade, down to 0.1 ppm as 
desirable in the bio-medical grade, The gamma photon activity as a function 
of time was calculated a3 detailed in Appendix C, The activity of the 236 Pu 
daughter nuclides was calculated from the growth of 1,91 year 220 Th. The 
calculation is detailed in Appendix D as a function of 230 Pu impurity, 
i.e, , 0,1 to 1.2 ppm. The gamma photons emitted by 23 6 Pu and the daughter 
nuelidoa arc given in Table 2-IV (1 ,9) . The rofaronne (9) data was uaod in 
this work. 

Alpha particles interacting with 18 0 may give rise to the formation of very 
short-lived compound nuclei which promptly decay accompanied by the 
emission of a neutron and gamma photons. The photon emission due to 

j 0 21 

o (a, n) Ne interactions is detailed in Appendix E and tabulated in 
Table 2-V. The radiated gamma photon fluxes produced by the * 8 0 (a, n) 2 *Ne 
reaction at time periods other than "fresh" were obtained by multiplying the 

fresh O («, n) Ne gamma yield values by the time decay factors given 
in Appendix E . 

It is indicated in reference (16) that the 17 0 abundance and (a, n ) cross 
section are approximately 10% of the corresponding values for 18 0. Though 
the O interaction has not been fully investigated and documented, and 
thus cannot be adequately treated at this time, it can be reasonably 
neglected. 

The gamma photon source spectra derived for each of the source categories 
I through IV above, are listed in Tables 2-II through 2-V. These sources are 
detailed for 0 (fresh), 1, 5, 10 and 18 years as a function of energy in 
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Tables 2-VI(a) and 2~VI(b) through 2-X(a) and 2-X(b). The (b) tables give 

236 

the sum of the isotopic photons and fission photons, and the Pu-chain 

18 

photons detailed in the (a) tables as well as the O reaction photons and 
their totals. Table 2 -XI summarizes these totals for each fuel age as a 
function of energy. 

238 (7) 

The calculations assumed the Pu half life as 87,4 years' 1 and the 

fission yield as 2.75 ± 0.01 neutron s/fiss ion ^ . The half-life of 241 Pu 

(7) 

was taken as 14,0 years . 

2 . 3 Gamma Photon Fluxes 

The gamma photon source spectra discussed in Section 2.2 were used to 
determine flux distributions at the RTG surface and exterior sRatial positions. 
The flux calculations were carried out for 1575, 3468 and 5679 w(th) model 
capsules. Figure 2 -II gives the dimensions of the three model capsules 
and their material zones as a function of thermal power. The material 
composition of the capsule zones are detailed in Table 2-XII 

Gamma photon flux transport to the RTG surface and exterior spatial positions 

(17) 

on the radial and axial midplanes were calculated using the ANISN and 

fl o\ (12) 

QAD ' computer codes. ANISN, a discrete ordinate code was used to 

(13) 

determine surface fluxes. QAD, an integrating point kernel code was 
used to calculate the fluxes at detector points outside the RTG. The 
gamma spectrum was assumed as unchanging from the capsule surface to 
detector positions distort from the capsule surface. The calculations were 
carried out for "fresh," 1, 5, 10, and 18 year old fuel. 
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The twenty group gamma photon cross sections were determined using a 

(19) 

modified version of the GAMLEG computer code. The cross sections 
in each of the energy intervals were averaged for PuO^ Product fuel with 
the gamma photon flux spectrum given in reference (1), 

Table 2 -XIII gives the calculated gamma photon surface flux as a function 

of energy and power, at the capsule radial midplane for fresh fuel. The 

fluxes are calculated for the decay of the plutonium Isotopes and daughters, 

236 

and fission but exclude Pu. The calculated values of k „ are 

effective 

noted as 0.17, 0.30 and 0.40 for the 1575, 3468 and 5679 w(th) capsules. 

IQ ^ 1 

Table 2-XIV presents similar photon flux data for the 0(a, n) Ne reaction 

as the photon source. Tables 2-XV through 2-XVII give similar photon flux 

236 

data as a function of fuel age for Pu and daughters as the photon source, 
for the 1575, 3468 and 5679 w(th) capsules, respectively. 

Tables 2-XVTII through 2-XXII give the calculated total photon fluxes at 
axial and radial detector positions at and beyond the capsule surface, 
i.e. , at cylindrical coordinates (r, 0 , z) as defined in Figure 2-III. The 
data in these tables corresponds to the total fluxes given in Tables 2 -XIII 
through 2-XVII at (r, 0 , z) = (capsule radial surface, 0, 0). The total 
fluxes given in Tables 2-XVIII through 2-XXII, Integrated over all sources, 
are shown graphically in Figures 2-IV through 2-IX as a function of detector 
distance r, with fuel age as the parameter. Figures 2-IV through 2-VI are 
for the radial midplane while Figures 2-VII through 2-IX are for the axial 
midplane . 

Figure 2-X gives the source-integrated total photon flux as a function of 
capsule power for r = 100 cm and 1 , 5 and 18 year old fuel. The curves 
demonstrate the ease with which datamaybe interpolated for other than the 
tabulated dependent variables. A further example is given in Figure 2-XI 
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where the total photon flux is plotted as a function of fuel age for the 
three model capsules. 

The calculations presented above allow the prediction of the gamma photon 

flux as a function of the following variables: 

236 18 

I fuel impurities such as Pu and Q, 

II fuel capsule power, and 

III fuel age. 

18 18 
If the O impurity is reduces from 100% to 10% as in the case of O 

18 21 

depleted oxygen then the O (a, n) Ne photon flux data may be multi- 
plied by the corresponding appropriate value in the range 1.0 to 0.1. 

2 36 

Again since the Pu and daughter photon flux data corresponds to a 

concentration of 1.2 ppm it may be readily corrected for other concentrations. 
236 

For example, if Pu is given as 0. 1 ppm then it is only necessary to 
236 

multiply the “ Pu fluxes" given in Table XV by the ratio 0. l/i .2 to 
get the correct photon flux. 

An example of the use of the calculated data presented in this report section 

is given in Appendix F. The example computes the photon flux at r = 100 cm 

2 3fi 

for a capsule power of 1575 w(th) , a Pu impurity of 1.6 ppm, "normal" 
oxygen, and 5 year old PuC^ fuel. 
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FIGURE 2-1 

TYPICAL PLUTONIUM FUELLED RTG CAPSULE 
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FIGURE 2 — II 

GEOMETRY OF MODEL FUEL CAPSULE 
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FIGUHt 2-V 

GAMMA PHOTON FLUX AS A FUNCTION OF Oi l I t : i OR RADIAL 
DISTANCE FROM 34tS8 w(tli) CAPStJLF AXIS WITH I" <WI : M mS i nr PARAMETER 
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F IGURE 2 VIII 

GAMMA PHOTON FLUX AS A FUNCTION * »l Of. 1 1 C'l'OM AXIAL 
| DISTANCE FROM 3460 with) CAPfiU I I M 1 1 -PI ANlJ 

WITH FUEL AGE AC Tilt PARAMETER 
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FIGURED IX 

GAMMA PHOTON FLUX AS A FUNCTION OF HI TFC I OH AXIAL 
DISTANCE FROM 6679 w(th) CAPSULE MIDP! ANF 
WITH FUEL AGE AS THE PARAMI.l EH 
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table z-n 

GAMMA PHOTONS FROM PECAY OF PLUTONIUM ISOTOPES 
ANP THEIR PECAY PRODUCTS 


U 1 

m luotopo 

Gamma Photon Cnorny 
— to-VJ- 

Ref . ( Q ) 

Abundance (»'of Uotcfdr; 

Baf-OD 

Hnfj.iMj 

• 1 

1 238 Pu 

17,0 

13,0< 3) 




43 .0 

3 , 8 x 1(T 2 

3,Hx I0~* 


*] 

1 

99,0 

8.0 x 1(T 3 

h.o x pr 1 

7.4 x 10" 3 

1 

1 

1 S2 . r > 

10-3 

1 . 1 X 1 o" 3 

0.0 x 1 0 “ 4 



207. H 

4.0 x 10 rb 

4,0 X |0“ f ' 

1 .2 x 1 0“ '' 



742 ,4 

9.0 x |(T 6 


4.3 x 10“'' 



709 .B 

3.9 X in-* 

9 , 0 X 1 0" 9 

2 . o x lo" ' 



7 HO , H 

0,0 x JO" 6 


2|H x 10“'' 

*) 


H 07 ,0 

0.0 x io” f) 

wiry '*fh'Ml 

0.0 x |0 1 

( 


891,3 



1.2 x In* 6 

' , . 


H79 , U 


2 .0 x 10“ ' 



HH2 . 9 



7,4 x m. 7 



920. » 



3,0 / 10 ' 

(■ 


04 1 . u 



...Ok .o 7 

’ t; 

i .• 


l oo | . 1 



0,3 x 10 * ' 


1 0-1 1 . ){ 
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17.0 
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/.Ox 10 i 
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3.0 X 10* 1 

1 . , X J0‘ ’ 


1 

■}. 

207.0 


4.4 x JO" 4 
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VI 0.0 


0.0 X 10* 4 



- 

17 3 . 0 

1 .2 x 10' 3 

1 . i X 10* * 


* 

/; 

41*1 .0 

1 .2 x 10* i 

11 , »J x 10“ 4 


.. /• i 


000.0 

H, Ox 10“ 3 



» B A 


770.0 

2.0 x 10"' 



• U I 
. ? 
: ■ 1 f 

: l 24 

17.0 
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4 >.3 
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090.0 

2 x 1 h ' 



•>,;f !• 

% t r 1 

34 V, 

1 00 . 0 


1 b 5 

1.21 x 1 o’ 4 * 

. V 

6 !" « 


14 9.0 

1 .0 X I0~ 4 

2.0 x 10" 4 

?■* 

ij • 24 '/a.-. 

00.0 

30.0 



t. 

! ,1 

101.0 

4.0 x 10~* 





200.0 

o X io“ 4 



IT 

Vi 9.0 

H X 1 0" 4 





37 0.0 

■1 X 10" 4 





003 . 0 

3 x 1 0" 4 



i 

1 

| 

772.0 

3 X 1 0~ 4 




I :u 

30.0 

14.0 
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r 

14 9.0 

1 .0 
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1 
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TABLE 2 -III 

GAMMA PHOTONS EMITTED FROM 
EQUILIBRIUM FISSION 


Gamma. Photon Energy 
(MeV) 

(15) 

Number of Photons/Fission 

Calculated 

Photons/Fission* 

0 . 1 — 0 . 4 

1.61 

1.61 

0.4-0. 9 

4.84 

4.86 

0.9-1.35 

0.50 

0.50 

1.35-1.80 

0.60 

0.61 

1.80-2.20 

0.31 

0.28 

2.20-2.60 

0.12 

0.11 

2.60-3.00 

0.01 

0.04 


* The photon yield per fission was calculated using analytic functions fitted 

to the above tabulated values of photon yield per fission. The functions 
used were: 


Gamma Ray Energy 
(MeV) 

. 1- .4 
.4-. 65 
.65-. 9 
.9-. 965 
,965-1.35 
1.35-1.60 
1.60-1.80 
1.80-3.0 


Function N (E) ** 
(photons/fission) 

21. 5E 
5.77e E 
22.60e~ 1,1E 
-121. 5E + 117.75 
0*38e 0,30E 
0.00045e 5,3E 
586.8e” 3,5E 
67.8e" 2,3E 
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TABLE 2"IV 

GAMMA PHOTONS PROM 236 Pu and DAUGHTER NUCLIDES 


Gamma Photon Energy 


Abundance (% of Isotope decay) 


236p u 

46, O^ 1 ' 

'(48.0)< 9) 4.7 X 10' 2 

3.1 x IQ" 2 


no*o 

1,2 x IQ" 2 

1.2 x 1U* 2 


165,0 

6,6 x IQ" 4 

6.6 x 10" 4 


520,0 


1.7 x 10" 4 


570,0 

, 

1.0 X IQ" 4 


645.0 


2.4 x IQ* 4 

212 Pb 

115.1 


0.7 


176.7 

82,0 (3) 

0.2 


238.6 

47.0 


300,1 


3.2 


415.2 


0.16 

2i2 Bi As a 

result of & decay (66% yield) 



40.0 


2.0 


288.0 


0.5 


460.0 


0.8 


727,0 

7.3 

7,1 


785,0 


l.i 


893.0 

6.6 

0.42 


953.0 

0,2 

0.10 


1074,0"] 

1079.0J 

' 0.66 

0.60 


1513.0 

0.86 

0.31 


1800.0] 

1809,0] 

small 

0.11 


1620.0 


1.8 

^* 2 Bi As a result of a decay (34% yield) 



288.2 


0.28 


328.0 


o.no 


434 .O' 




453.0 


0.42 


473.0 




493 « 0^ 



208 T £ 

280.0 

10.0 



511.0 

25.0 

23.0 


583.0 

80.0 

86.0 


860.0 

15.0 

12.0 


2614.0 

100.0 

100.0 


The gamma activity ox 236 Pu WQg calculated as discussed in Appendix C, 

The radioactive decay chains and corresponding daughter nuclide activity are 
calculated in Appendix D . 


♦From reference (3) 
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TABLE 2-V 

GAMMA PHOTONS DUE TO l8 0(a ,n) 2l Ne REACTION 


Gamma Photon Energy* 


(MeV) 


0.35 

4 x 10 3 

1.38 

8.9 xlO 2 

1.90 

1.8 x IQ 2 

2.40 

1.8 x 10 2 

2.70 

1.8 x 10 2 


♦Although references (10) and (11) suggest 0.35, 1 .38, 1.75 and 2.87 gamma 
photons, it is felt that the 1.75 MeV and 2.87 gamma rays contained the 
above as unresolved resonances . 

The photons emitted per neutron emission-, neutron emission rates and other 
pertinent data are presented and calculated in Appendix B. 
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Includes prompt and fission product gamma photons 
1.2 ppm 



GAMMA PHOION ACTIVITY AT TIME TEN YEARS (T = 3 0) 
(Photons/'sec-gm PuC> 2 ) 





* Includes prompt and fission product gamma photons . 
** 1 .2 ppm 
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TABLE 2 -XII 


MATERIAL PROPERTIES OP 
MODEL FUEL CAPSULE 


Density 

Zone Material gm. /cm 3 

1 Void . 0.0 

2 Tungsten (W) 19.3 

3 Plutonium-Fuel(Pu 02 ) 7.306*** 

Plutonium (Pu) 

Oxygen (O) 

4 Tantalum(Ta) 16.6 


Atomic Density* 
atoms/cm 3 

Modified 
Atomic Density 

0.0 


.063190 

4.67628 

0.01632 

1.79520 

0.01632 

1 .53408 

0.03264 

0.26112 

0.05526 

4.03380 


* This density corresponds to the P Q Legendre expansion term in the Klein-Nishina 
cross sections. 

** This density corresponds to the Pj_-P 5 Legendre expansion terms in the Klein-Nishina 
cross sections. The P 1 -P 5 Legendre terms of hydrogen were utilized and corrected 
by the Z (atomic number) of the isotope used in the analysis . 

*** The Plutonium Dioxide fuel density, 10.7 gm /cm 3 , was multiplied by the volume 
fraction 0 . 683 . 



Total 1 .95648E 05 2.65065E+05 3.D4190E+05 
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GAMMA PHOTON FLUX (Photon/cra 2 -sec) FOR 1575 w(th) CAPSULE AS 
A FUNCTION OF DETECTOR POSITION AT NOTED FUEL AGE; 
SOURCE — 236 Pu AND DAUGHTER^ 
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TABLE 2 -XXII 

GAMMA PHOTON FLUX (Photons/cm 2 -sec) FOR 5679 w(th) CAPSULE AS 
A FUNCTION OF DETECTOR POSITION AT NOTED FUEL AGE; 
SOURCE — 236 Pu AND DAUGHTERS 
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3. 


NEUTRON RADIATION CHARACTERISTICS 


3.1 Introduction 


This report section presents the result of calculations of the neutron flux 
emitted by plutonium dioxide fuelled RTG capsules. Neutron flux dis- 
tributions were determined as a function of energy for the three RTG 
model capsules already described in Section 2, namely: 1575, 3468 
and 5679 w(th) , Thaflux calculations were based on neutron yield 
spectra which were derived to account-for plutonium spontaneous fission 
and both oxygen and contaminant light element (a, n) reactions. 

Since the RTG .emitted neutron flux is a function of the total neutron yield 

per gram of give a plutonium dioxide , _the yield spectra must be obtained 

as discussed in Section 3.2. Measurements and calculations of the 

spectrum of alpha particles and of neutrons in production plutonium 

(1-3 >2-27) 

dioxide fuel are presented, in the literature ' . Unfortunately, 

there is considerable variation, with respect to the emitted spectral 
components and their intensities. 

The total neutron yield for plutonium dioxide with oxygen of natural 

abundance and no impurities or self multiplication may be as low as 
4 238 (28^ 

1.7 x 10 neutrons/sec-gm of Pu . Plutonium dioxide sources 
larger than 100 w(th) , typically have a total specific neutron yield between 

A A n on /i o q\ 

2.0 x 10^ and 3.0 x 10* neutrons/sec-gram Pu ' . The total 

18 21 

specific yield may be taken as 75% due to the 0(<*, n) Ne reaction, 

238 

12% due to the spontaneous fission of Pu and 10% due to the con- 
taminant light element (fit, n) reactions with the remaining 3% due to lesser 

17 21 

known phenomena such as 0(<*, n) Ne, (y, n) , etc. Typical yield values 


I. Plutonium fission — — spontaneous fission neutron yields 
of-2586 (± 398) neutrons/sec-gm 238 Pu have been reported^ 9 ' 3 °^ , 

II. Oxygen { a , n ) reactions 17 0/ 18 C (a, n) 21 Ne neutron 

yields ranging from 12, .400 to 14,500 neutrons/sec-gm Z38 Pu., 
have boen reported in references (1) and (28), respectively, 

III. Impurity (a, n) reactions an impurity (a, n) neutron . 

yield of 4000 neutrons/sec-gm Pu is typical for commercial 
fuel, though values as high as 10,000 neutrons/sec-gm 238 Pu 
have been reported. 

There are current efforts ^ to produce PuC> 2 depleted in 17 0 and 18 0 in 
order to reduce the neutron emission rate to about 20% that of the natural 
oxygen .PuO^ . 

Section 3.2.1 describes the method used to determine the neutron fission 

spectrum. Section 3.2.2 and Appendices G, H, I and J detail the method 

18 2.1 

used to determine the 0(a, n) Ne neutron spectrum. The impurity 
element (a, n) spectrum was calculated as described in Section 3.2.3 
and Appendix K. 

The neutron yield spectra were used as source terms for determination of 

neutron transport to the fuel capsule radial and axial surfaces, and surface 

flux spectra thus obtained. Neutron transport was calculated by the method 

( 12 ) 

of discrete, ordinates . Neutron flux geometric reduction from the capsule 

surface to exterior spatial locations was calculated by means of the removal 

(13) 

cross section technique and point kernel geometry . These calculations 
and obtained results are discussed in Section 3.3. 
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3.2 Neutron Souroe Spectra 


3,2.1 Fission Neutrons 

The neutron fission spectrum was determined from the reported neutron 

total spontaneous fission yield rate of 2.6 x 10 neutrons/seo-gm Pu' . 

The 2.6 x 10 3 value does not include the spontaneous fission of 240 Pu 

because of its low abundance in the fuel and the fact that its fission rate, 

29 neutrons/ sec-gm 2 ^Pu, is less than the uncertainties indicated in the 

238 

literature for the spontaneous fission rate of Pu. The fission. spectrum 

(15) 

was. determined from the uranium fission spectrum relationship normalized 
to the plutonium total spontaneous fission yield, as 

N(E) = 2.00154 x 10 3 E„ 1/,Z .exp (-E /l .29) , (3.1) 

S S3 

where E is the fission neutron energy. An alternative relationship has 
s 

been reported which gives the number of spontaneous fission neutrons 

238 ( 241 

produced per MeV per gm of Pu as' ^ 

N(E ) = 2.04 x 10 3 .E 2 . exp (-E /1.34) (3.2) 

s s s 

This relationship is based on a total spontaneous fission yield rate of 
3 238 (28) 

2.8 x 10 neutron/s ec-gm 4 Pu ' . The total number of spontaneous 

238 

fission neutrons produced by Pu in any energy interval may be obtained 
after integration of equation (3.1) , over the energy interval, as 

r% f E tt 

N = I N(E s )dE s = 2.00154 x 10 3 I E g 1//2 exp (-E s /1 .29) dE s , (3 .3) 

where, E L and E H are the lower and upper limits of the neutron energy interval 
in energy units of MeV . 
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I 

I 


The total spontaneous fission neutron yield rate was calculated for 

twenty-three (23) energy intervals, in the range 0.025 eV to 10 MeV. 

The results of the calculation are given in Table 3-1 for the RTG plutonium 

dioxide fuel. In Table 3-1 column 2 is the spontaneous fission yield- per 
2 3 B 

gm of The spontaneous fission yield in each energy interval was 

divided by the total neutron fission yield to determine the spontaneous 
fission spectrum given in column 5 of Table_3=I, . 


In the neutron flux results presented later in this section, the value of 
^eff corres P° n< ^ s the ratio of the calculated total neutron fission rate 
in the subcritical system volume to the assumed neutron fission source 
normalized to unity in accord with the group distribution given in Table 
3-1. The fission rate for each neutron group g, and for each volume 
element, was integrated as: 

. (3.4) 

2 th 

where N. is the scaler flux (neutrens/cm -sec); V. is the i elemental 
volume (cm ); and vxr. is the number-of neutrons released per fission 
multiplied by the fission cross-seotion corresponding to group g. A 

more extensive multigroup treatment of-the multiplication constant, k 

(17 31) ett 

etc. is given in the references ' 

3.2.2 Oxygen («, n) Neutrons 

The oxide in the plutonium dioxide fuel source contains natural oxygen 

16 17 18 

which consists of the oxygen isotopes O, O and O in the relative 

( 23 ) 

abundances 99.759, 0.037 and 0,204%' . Alpha particles interacting 
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with low atomic number isotopes may give rise to the formation of very 
short-lived compound nuclei which, promptly decay accompanied by the 

emission of a neutron. The aross section for this (ft, n) reaction is negligible 

16 17 90 1ft 

for the relatively stable O isotope, However 0(ft, n) 4 Ne and * Q(ft, n) 

21 

Ne reactions represent the principal sources of neutrons in plutonium 
dioxide . 

lft 21 

The total neutron yield produced by the 0(ft, n) Ne reaction in the 

plutonium dioxide product fuel are tabulated in the current literature 
3 4 21 32 — 35 ) 

' ' ' ' ' . The (a, n) neutron yields in the references are 

expressed as the numerical difference between the measured total neutron 
yield and the calculated neutron fission yield. Variations in the (ft, n) 
spectrum result from uncertainties in the calculated fission yield and 
lack of knowledge of light element contaminations in the plutonium dioxide 
fuel. Though oxygen (a, n) neutron yields have been reported^ 1 ' 24 ^ 
information has not been given in sufficient detail to allow spectral 
calculations to be carried out. In the remainder of this subsection and 
Appendices G, H, I and J the (a, n) total neutron spectrum and the methods of 
calculation are presented. 

In the present work, the (a, n) neutron yield from the irradiated oxygen 

1 8 

nuclei was assumed to be totally due to the O isotope since the maximum 

17 20 18 21 
yield for the 0(ft, n) Ne reaction is ~ 10% of the O(ot, n) Ne reaction 

yield^ 1 ^ the abundance is 0.037% as compared to 0.204% for * 8 0. 

17 21 

Further consideration of the 0(ft, n) Ne neutron yield is not possible 
until more experimental information becomes available. 

An alpha particle of initial energy E, collides with the atoms of the trans- 
port medium and is thus degraded in energy until thermal equilibrium with 
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the medium is achieved. Each alpha particle collision with the 18 0 
nucleus correspond!:, to an inelastic (a, n) reaction providing the alpha 
energy is greater than the threshold energy. The threshold is the lowest 
energy for which the (a, n) reaction can occur. The alpha particle "slowing 
down" process gives rise to an energy spectrum of inelastic collision 
neutrons. If isotropic scattering in the center of mass coordinate system 
is assumed, tho number of noutrons with energy E n , generated by tho alpha 

1 Q 

particle collision with Q target nuclei may be obtained as 

N(E n ) dE n - N (6 crn ) ein^ m d8 cm 
2 

where, ® cj ^is the scatter angle in the center of mass coordinate system. 

The reader may refer to texts for a discussion on the penetration of charged 
particle si^. 

22 i_o 

The intermediate Ne nucleus resulting from the 0(«, n) reaction decays 
to various ^Ne levels as follows (2 1 ' 32 ' 33 

55% to the ground state, 

35% to the first excited 0,35 MeV state, 

^ 10% to the second excited 1.73 MeV state, 

~ 1% to the thir’ excited 2.84 MeV state. 

The (a, n) reactions were calculated for each of these inelastic collision 

schemes. The Q values for the reactions were taken as -0.70, -1.05, 

-2.43 and -3.54 MeV for ground, first, second and third states, respectively. 

238 ( 1 

The two alpha particles emitted by Pu , namely; 


Energy (MeV) 
5.495 
5.432 


Alpha Particle per Disintegration 
0.72 


0.28 

were taken as a single 5.50 MeV alpha particle for all calculations. 

The dynamics of the (a, n) 21 Ne reaction are given in Appendix G. 

The alpha particle energy loss per unit path length, (dEe/dx)*, resulting 

from collisions with oxygen and plutonium atoms 4n the fuel, is discussed 

in Appendix H. The G (tt, n) No neutron yield spectrum is discussed 

in Appendix I. Appendix J describes the numerical integration of tho * 8 0 
2 1 

(a, n) Ne spectrum over the twenty-throe energy groups used in tho 
present work. 

18 21 

The O (a, n) Ne neutron spectrum was determined for the PuO product 

fuel in twenty-three energy groups in the range 0.025 eV to 10.0 MeV. 

Table 3-II presents the results of the calculations. The data are normalized 

to a neutron yield Value of 1.22 x 10 4 neutrons/sec-gm 238 Pu; an experimentally 

measured value of 1.24 x 10 4 neutrons/sec-gm 238 Pu has been obtained 
17 1R 

for reactions with O O w ". Table 3 -III is the spectrum normalized 

to unit emission. 

3.2.3 Impurity (Ot, n) Neutron Yield 

Commercial grade PuO^ contains small quantities of uniformly mixed low 
atomic number impurity elements which yield neutrons through the (a, n) 
reaction. The neutron yield# which may be significant, varies with each 
commercial fuel feedstock. Thus in order to calculate the total PuOg 
neutron yield it is necessary to explicitly identify the impurity elements 
present. Since the yield spectra have not been generally reported they 
were assumed, in the present work, to be similar to the 18 0 (a, n) 21 Ne 
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yield spectra as given in Table 3-III. This assumption if the impurity 
concentration is relatively small and better data is unavailable. 

Table 3-IV presents an analysis of a typical Pu(N0 3 ) 4 feed stock used in 
the production of commercial PuO^ fuel^ 38 \ The 338 Pu isntope emits 
~5,5 MeV alpha particles which in turn interact with the impurities. The 
5.5 MeV alpha particles have sufficient energy to exoeed the (a, n) 
threshold and the coulomb barrier of the low atomic numbered impurities . 
Neutron yield data for the major impurity elements are determined in 
Appendix K; the yield units ares neutrons/sec-^jm Pu per ppm of impurity. 

For example if a Li impurity of 5 ppm (by weight) is present in Pu0 2 , 
then the neutron yield will range from,, 5 x 4.3 = 21.5 to 5 x 5.5 = 27.5, 
neutrons/sec-gm Pu. The range variation results from the range distribu- 
tion of experimental data. An excellent compilation of light element 
(a, n) reactions is given in Chapter I of reference (40); reference (41) 
is also recommended. 

3.3 Neutron Fluxes 


The neutron source spectra yields discussed in Section 3.2 were used to 
determine neutron flux distributions at the RTG surface and exterior spatial 
positions. The flux calculations were carried out for three fuel capsules 
already described in Section 2, namely: 1575, 3468 and 5679 w(th). 

Figure 2 -II gives the dimensions of the three model capsules and the material 
sones as a function of power. The material composition of the capsules' 
zones are detailed in Table ?. -XIII. 
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The neutron transport to the RTG surface and exterior spatial positions 

(17) (18) 

was calculated using the ANISN' ; and QAD' ' computer codes. 

ANISN was used to determine surface fluxes and QAD was used to obtain 
exterior fluxes. The neutron collision group constants were calculated 

( 42 ) * 

using the "Evaluated Nuclear Data File"' ' and the"R, J. Howerton 

(43) 

Evaluated Data File" . The microscopic collision cross sections were 

averaged for twenty-three neutron groups with a neutron flux spectrum 

thermalizsd to 20°C and continuously varied as a function of energy E, 

238 ( 1 ) 

up to approximately 0.2 MeV. The measured PuOg' ' ern * 8S * on neutron 

spectrum was assumed as the fission weighted flux spectrum up to 7 MeV. 

238 

From 7 MeV to 10 MeV, + he Pu fission spectrum was assumed for 
weighting. Additional calculations showed that even when thermalized 
up to 600°C the thermal neutron fluxes remain relatively unchanged. 

Table 3-V gives the calculated neutron flux as a function of energy and 

18 21 

power at the capsule radial midplane surface for O (tt, n) Ne and 

fission neutrons. Table 3-VT gives similar radial surface fluxes for fission 

neutrons only. The calculated values of k .. are noted as 0.17, 

effective 

0. 30 and 0,40 for the 1575, 3468 and 5679 w(th) capsules. 

Tables 3 -VII, 3 -VIII and 3-DC give the calculated total neutron flux at 
axial and radial detector positions at and beyond the capsule surface, 

1. e. , at cylindrical coordinate positions (r,8,z) as defined in Figure 2-III. 
The data in these tables corresponds to the total fluxes given in Tables 
3-V and 3 -VI at (r,8,z) = (capsule radial surface, 0,0). The total 
fluxes given in Tables 3-VII to 3-K are shown graphically in Figures 3-1 
through 3-IV as a function of detector distance r, with capsule power as the 
parameter. Figures 3-1 and 3-II are for the radial midplane while 

Figures 3-III and 3-IV are for the axial midplane. Figures 3-1 and 3— III 
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18 21 

are for O (a, n) Ne plus fission neutrons while Figures 3-II and 3-IV 

18 9 7 

are fission neutrons only. Data for O (a, n) Ne may be obtained by 
subtraction since the values of k are less than unity, i.e. , since the 
fuel capsules are subcritical. 

The surface neutron flux produced by the impurity light element (a, n) 
neutrons may be determined as 



where 

* o = surface neutron flux produced by the oxygen (a, n) reaction, 
eg. Table 3-V minus Table 3-VI; 

Y 0 = neutron yield produced by the oxygen (a, n) reaction per 

second per gram of Pu, eg. 9882 n/sec-gm Pu based on 80% 
238 Pu in Pu and 1.22 x 10^ n/sec-gm 338 Pu; 

Yj » neutron yield produced by the impurity element (<*, n) reaction 
per second per gram of Pu, eg. the Appendix Table K-II data 
gives Y^ = 191 x 10 = 1910 n/sec-gm Pu for 10 ppm of boron 
(using the maximum yield per ppm of impurity value) . 

Equation (3.5) assumes that the light element (a, n) neutron spectrum 

18 21 

corresponds to the O (a, n) Ne neutron spectrum. 

An example of the use of the calculated data presented in this report 
section is given in Appendix F. The example computes the neutron flux at 
r = 100 cm for a capsule power of 1575 w(th) , "normal" oxygen and a 
100 ppm boron impurity. 
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FIGURE 3-1 

NEUTRON FLUX AS A FUNCTION OF DETECTOR RADIAL DISTANCE 
FROM CAPSULE AXIS WITH POWER AS THE PARAMETER 
(NEUTRON SOURCE: 18 0(«,n^ Ne AND FISSION) 
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FIOURK 3 -11 

NEUTRON FLUX AS A FUNCTION OF DETECTOR RADIAL DISTANCE 
FROM CAPSULE AXIS WITH POWER AS THE PARAMETER 
(NEUTRON SOURCE: FISSION) 
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FLUX AT DETECTOR (NEUTRQNS/CM Z -SEC) 



DETECTOR AXIAL DISTANCE (CM ) 





FLUX AT DETECTOR (NEUTRONS/CAT-SEC) 


1 


FIGURE A IV 

NEUTRON FLUX AS A FUNCTION OF DETECTOR AXIAL DISTANCE 
FROM CAPSULE AXIS WITH POWER AS THE PARAMETER 
(NEUTRON SOURCE: FISSION) 
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TABLE 3-1 


PLUTONIUM ISOTOPE SPONTANEOUS FISSION NEUTRON SPECTRUM 


Neutron Energy 
. Interval 

Neutrons 

(MeV) 

sec-gm^8p u 

0.025 - 0.10 x 10" 6 

3.72 x IQ" 8 

0.010 - 0.30 x 10“ 5 

6.90 x IQ" 6 

■0.030 - 0.30 x 10" 4 

2.13 x 10" 4 

0.030 - 0.555 x 10" 3 

1.73 x 10~ 2 

0.555 x 10" 3 - 0.0170 

2.92 

0.0170 - 0.0449 

9.52 

0.0449 - 0.122 

41.39 

0.122 - 0.201 

55.97 

0.201 - 0.331 

108.94 

0.331 - 0.546 

202.24 

0.546 - 0.702 

152.14 

0.702 - 0.90 ' 

190.74 

0.900 - 1.16 

237.86 

1.16 - 1.49 

272.56 

1 -49 r 1.91 

294.11 

1-.91 -2.45 

295.72 

2.45 - 3.14 

266.52 

3.14 -4.04 

214.13 

4.04 - 4.46 

64.58 

4.46 - 5.18 

76.26 

5.18 - 6.66 

76.51 

6.66 •• 8.55 

30.93 

8.55 - 10.0; 

7.00 

Total 

2.6 x 10 3 


JIM U U U UUD 

sec-gm PuC >2 

Neutrons/seo-cm 

(PuC 2 Fuel) 

3 Fission 
Spectrum 

2 . 64 x 

10" 8 

19.3 x 10** 8 

0.0 

4.93 x 

10" 6 

36.03 x IQ" 6 

0.0 

1 . 52 x 

Hf 4 

11.11 x 10" 4 

0.0 

1 .24 x 

IQ" 2 

9.06 x IQ -2 

• 0.0 

2.09 


15.27 

.00113 

6.80 


49.70 

.00368 

29.55 


216.0 

• .01598 

40.0 


292.32 

.02162 

77.78 


568.43 

.04204 

144.39 


1055.2 

.07805 

108.62 


793.8 

.05871 

136.18 


995.22 

.07361 

169.82 


1241.07 

.09180 

194.59 


1422.08 

.10518 

210.0 


1534.7 

.11351 

211.12 


1542.9 

.11412 

190.28 


1390.59 

.10286 

152.88 


111 7 . 27 

.08264 

46.11 


337.0 

.02493 

54.44 


397.86 

.02943 

54.62 


399.17 

.02952 

15.77 


• 115.25 

.00853 

5.00 


36.54 

.00270 

1.86 x lO' 3 

13.52 xlO 3 

1.0000 
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TABLE 3-II 

TOTAL NEUTRON EMISSION FROM THE 1 8Q(a,n) 21 Ne 
REACTION IN Pu0 2 FUEL* 
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volume fraction 0.693 was used for this study 


TABLE 3 -III 


18 0(a,n) 2l Ne NEUTRON SPECTRUM 


Neutron Energy 
Interval 
(MeV) 


Noutron Spectrum* 
(N(E) 


2.5 x 10~ 8 ~1Q’7 

0.0 

10” 7 -3 . 0 x IQ’ 6 

0.0 

3.Q x lQ’ 6 -3.0 x 1Q“ S 

0.0 

3.0 X 11T 5 -S.5S x IQ’ 4 

0.0 

5.55 x 10~ 4 - 1,70 x 10’ 2 

.00022 

1.70 x 10~ 2 -4.49 x 10’ 2 

.00045 

4.49 x 10’ 2 -1 .22 x 10’ 1 

.0020 

1.22 x 10 _1 -2 ,01 x 10" 1 

.00216 

2.01 x 10" l -3.31 x 10“* 

.00376 

3.31 x 10” 1 -5 .46 x 10- 1 

.0078 

5.46 x 10*" 1 -7 . 02 x 10’ 1 

.00821 

7.02 x lO’^S.QO x 10" 1 

.01019 

9.00 x 10~ 1 -1 . 16 

.02883 

1.16 - 1.49 

.04875 

1.49 - 1.91 

.07483 

1.91 - 2.45 

.17095 

2.45 - 3.14 

.34745 

3.14 - 4.04 

.26938 

4.04- 4 .46 

.02475 

4.46 - 5.18 

.000105 

5.18 - 6.66 

.0 

6.66 - 8.55 

.0 

8.55 - 10.00 

.0 

Total 

* Number of neutrons emitted 
per unit neutron emitted by 
18 0(a , n) 21 Ne reaction 

0.99983 
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NEUTRON FLUX AS A FUNCTION OF DETECTOR POSITION 
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4. 


DISCUSSION OF RESULTS 


Basic gamma photon and neutron source charactistics have been determined 
and flux data obtained for plutonium dioxide RTG fuel capsules. The fuel 
capsules were hypothesized based on the SNAP-27. The methods for 
characterizing the sources have been given in adequate detail. The cal- 
culation of the photon and neutron fluxes were carried out using computer 
codes in the public domain. The source data may be used as input to 
transport calculations in other than SNAP-27 configurations, as for example 
in SNAP-19 where the fuel is not annular but consists of solid discs, 
i.e. , "hockey-pucks. " Alternately the flux data in this report may be 
corrected to approximate data for other fuel configurations. 

The results may be considered as useful starting points for analyzing RTG 
radiation interference With science experiments or instrumentation on 
space probes. Similarly it may be used to determine biological hazards 
in the presence of RTG’s. The results may be modified as indicated for 
analyses of biomedical problems where PuOg power sources are proposed, 
as for example in "heart programs." 

Although the data in Sections 2 and 3 do not include exposure dose information 

this may be readily obtained by employing standard conversion factor 
(44 45) 

information ' . Table 4-1 presents RTG capsule radial gamma photon total 

dose rate at one meter as a function of fuel age with power as a parameter. 

The table also indicates the percentage contribution from the source com- 
2 36 1ft 

ponents: Pu, O, and Pu isotopes and daughters. It is noted that the 

isotopes and daughters are the major source (~ 97%) component in fresh 
fuel, whereas Pu is the major source in aged fuel (~ 91%). Figure 4-1 
gives the one meter radial and axial dose rate as a function of fuel age for 


76 


the 1575 w(th) model capsule; this size approximates SNAP-27. Reference 
(46) gives radial and axial experimental values of 8 and 1,6 mr/hour for 
a two year old 1500 w(th) SNAP-27. Allowing for the 5% difference in source 
power, this value is in very good agreement with the present calculations 
as seen in Figure 4-1. 

The capsule surface neutron flux data given in Table 3-V and the reference 

(45) conversion factors were usad-to-obtain the surface dose rate for a 

1575 w(th) capsule. The data in Table 3 -VII was used to convert surface 

dose rates to one-meter radial and axial dose rates of 52.4 and 30.5 mrem/hr. 

Since these dose rates are based on source yield data which accounts only 

for O (a, n) Ne reactions and fission,. ..a '.correction' for (Of, n) impurity 

reactions is required. If impurities concentrations are not identified but * 

total yield data is , an estimation of the (Of, n) impurity flux may be carried 

out. Reference (24) gives total neutron yield for SNAP-27-1 as 2.2 ± 0.1 x 
4 238 

10 neutrons/sec-gm Pu, and the dose rates above (52.4 and 30.5 mrem/hr) 

4 238 

are based on 1.48 x 10 neutron/sec-gm Pu, which suggests an impurity 
4 238 

yield of 0*72 x 1£) — rxeutxon/sec-gm Pu, ie. , an additional ~ 49% yield. 

For this yield radial and axial total neutron dose rates of 78 and 45 mrem/hour 

are predicted. The actual total neutron yield value depends on the PuO^ 

quality control during the production process. If the total neutron yield 

238 

is taken as 19,000 neutrons/sec-gm Pu, the SRP value as in reference 
(37), then radial and axial total neutron dose rates of 64 and 38 mrem/hour 
are predicted. These dose rates may be compared with references (47,48) 
radial and axial measurements of 93 and 42 mrem/hour, respectively. The 
reference data represents a mean for three different capsules. 

Since the experimental dose data of the references includes laboratory 
neutron and photon scattering it is suggested that this may be the major 


reason for the difference between experiment and prodiction. The backscatter 
dose is a function ot experiment geometry, structural material composition 
and thickness as well as neutron and photon spectral distribution. For 
gamma photons the increase in excess of the direct dose, resulting from 

( 29 ) 

concrete floor and wall scatter, would be only ~ 5% . It is possibly 

for this reason that the above gamma photon predictions are in good agree- 
ment with experiment. The increase in excess of the direct dose for 
neutrons is considerably greater for concrete floor and wall scatter, with 
a value of " 30% estimated as probable; this percentage is based on 


estimates using references (50) and (51) albedo data. In order to obtain 
an accurate determination of backscatter the details of the specific 
experiment are required. A '30% backscatter' for neutrons would increase 
the radial and axial total direct-dose rates of 64 and 38 mrem/hours to 
yield 'laboratory dose' rates of 83 and 49 mrem/hour respectively. These 
dose rates are within the experimental data accuracy of 93 and 42 mrem/hour 


± 20% 


(47) 
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APPENDIX A 

GAMMA ACTIVITY DUE TO DIRECT DECAY OF 
241 

Pu AND DAUGHTER NUCLIDES 

Ths direct decay of 241py yields a 145 koV gamma photon whose intensity 
is obtained as: 

photons/sec ™gm Pu0 2 = 3.92x10*2 dls x ( . 008) x ( . 88 1) x (abundance) 

soc-gm 241p u 

where 


and 


0.008 is the fraction gm 241p u 

gm Pu 


0.881 is the fraction gm Pu 

gm Pu 02 

Thus, the 145 keV photon activity is obtained as 


* 3.92x10*2 dis x ,008x .881x1.6 x 10" 6 
sec-gm ^ 4i Pu 

= 4.42 x 10 4 photons (145 keV) 
sec-gm Pu 02 


There are two decay modes available to 241p u# The gamma photon yields 
241 

of the Pu daughter isotopes in these chains is presented in the remainder 
of this appendix. The decay chains available to 241p u are: 


241 


Pu 


JL* 

14y 


241 


Am 


a 


237 


432y 


N 


P 2 . 14 x 10 6 y 


(I) 


241 


Pu 




J_ 

6.75day 


■* 



2 . 14 x 10 



. (II) 


Decay Chain 1 

Decay chain I gives rise to a gamma spectrum which is calculated from the 
following abundances 


Gamma Energy (koVl 


Abundance (% of isotope decay) 


60 

36. 

101 

4 x 10" 2 

208 

6 x 10" 4 

335 

8 x 10" 4 

370 

4 x 10” 4 

663 

5 x 10 -4 

772 

3 x 10‘ 4 


The gamma activity may be calculated from the Bateman equation. Where 
the decay chain is represented as: 


A 



B 



* 


the activity of 24 ^-Am, 


represented by BX q , 


is obtained as 


I 


where 

= dis/sec, of daughter nuclide B, 

Naq ® initial number of atoms A; equal to Avagadros number 

(N 0 ) divided by the gm-atomic mass of the parent A, given 
as M a 


The decay constants Xa • ^b oorres P ondin 9 to to 0 parent and daughter 
nuclides were calculated os, ^ «= .0495 year = 1.567 x 10 sec 


^ = .001604 year"'* 1 


24 1 

The activity of Am is thus obtained as 


dis 

f qm x 

f No\ 

Xa Xg \ 

sec-gm Pu0 2 ! 

^gm Pu0 2 I 

I m a | 

iwl 


(e"^At_ e "^Bt) x (abundance) 


Decay Chain II 


Decay...chain II yields four photon lines which may be calculated from the 
relationship; 


photons 

= (Xa No 


Um 241 Pu\ 

sec-gm Pu0 2 

m A 

e“Vtj 

^gm Pu0 2 ^ 


x (abundance) x (yield) 


where the yield is assumed as 2.3 x 10“ 5 
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APPENDIX B 

TABULATION OF PROMPT FISSION AND 
EQUILIBRIUM FISSION PRODUCT GAMMA PHOTONS 


This appendix presents tabulations of prompt fission and equilibrium fission 
product gamma photon emission. The tabulations are given in support of 
Section 2 of this report . 


I 


TABLE B-I 

PROMPT FISSION AND EQUILIBRIUM FISSION 
PRODUCT GAMMA PHOTONS* 


Gamma Ray Energy 
Interval (MeV) 

7 . 0 - 6.0 

6 . 0 - 5 . 0 . . 

5 . 0 - 4.0 

4 . 0 - 3.0 

3 . 0 - 2.0 

2 . 0 - 1.8 
1 . 8 - 1. 6 
1 . 6 - 1 . 4 
1 . 4 - 1. 2 
1 . 2 - 1.0 
1 . 0 - 0. 9 
0 . 9 - 0. 8 
0 . 8 - 0. 7 
0 . 7 - 0. 6 
0 . 6 - 0. 5 
0 . 5 - 0.4 
0 . 4 - 0. 3 
0 . 3 - 0. 2 
0 . 2 - 0.044 

0 . 044 - 0.001 


Prompt Fission 
Gamma Spectrum 

0.007 

0.020 

0.060 

0.179 

0.538 

0.198 

0.247 

0.308 

0.384 

0.478 

0.302 

0.380 

0.479 

0.602 

0.758 

0.954 

1.201 

0.0 

0.0 

0.0 


Equil, Fission Product 
Gamma Spectrum 

0.0 

0.0 

0.003 

0.027 

0.267 

0.173 

0.312 

0.267 

0.117 

0.106 

0.307 

0.888 

0.9910 

1.077 

1.000 

0.905 

0.753 

0.538 

0.323 

0.0 


Total Fission 
Gamma Spectrum 

0.007 

0.020 

0.063 

0.207 

0.815 

0.371 

0.559 

0.575 

0.501 

0.584 

0.609 

1.268 

1.470 

1.679 

1.758 

1.859 

1.954 

0.538 

0.323 

0.0 


* Tabulation units are: "Number of photons emitted/fission." 


TABLE B -II 


PROMPT FISSION AND EQUILIBRIUM FISSION 
PRODUCT GAMMA PHOTONS 

Total Fission Photons 


Gamma Ray Energy 
Interval (MaV) 

Gamma Spectrum 
(y /fission) 

sec 

keff=0 ,17* 

-gm PuC >2 
0.30* 

0.41* 

7. 0-6.0 

0.007 

5.0 

10.0 

15.1 

6. 0-5.0 

0.020 

14.0 

28.2 

43.0 

5. 0-4.0 

0.063 

45.0 

88.8 

1.355 x 10 

4. 0-3.0 

0.207 

147.0 

291.9 

4.451 x 10 

3. 0-2.0 

0.815 

5.787 x 10 2 

11.492 x 10 2 

1.752 x 10 

2. 0-1. 8 

0.371 

2.634 x 10 2 

5.231 x 10 2 

7.977 x 10 

1.8-1. 6 

0.559 

3.969 x 10 2 

7.882 x 10 2 

1.202 x 10 

1.6-1. 4 

0.575 

4.083 x 10 2 

8.108 x 10 2 

1.236 x 10 

1.4-1. 2 

0.501 

3.557 x 10 2 

7.064 x 10 2 

1.077 x 10 

1. 2-1.0 

0.584 

4.147 x 10 2 

8.234 x 10 2 

1.256 x 10 

1.0-0. 9 

0.609 

4.324 x 10 2 

8.587 x 10 2 

1.309 x 10 

0.9-0. 8 

1.268 

9.003 x 10 2 

17.879 x 10 2 

2.726 x 10 

0.8-0. 7 

1.470 

10.044 x 10 2 

20.727 x 10 2 

3.161 x 10 

0.7-0. 6 

1.679 

11.921 x 10 2 

23.674 x IQ 2 

3.610 x 10 

0.6-0. 5 

1.758 

12.482 x 10 2 

24.788 x 10 2 

3.780 x 10 

0.5-0. 4 

1.859 

13.200 x 10 2 

26.212 x 10 2 

3.997 x 10 

0.4-0. 3 

1.954 

13.874 x 10 2 

27.552 x 10 2 

4.202 x 10 

0.3-0. 2 

0.538 

3.820 x 10 2 

7.586 x 10 2 

1.157 x 10 

0.2-0.044 

0.323 

2.294 x 10 2 

4.554 x 10 2 

.695 x 10 

0.044-0.001 

0.0 

0.0 

0.0 

0.0 


* The fission rates determined for this study^^ were 710, 1410 and 2150 
n/sec-gm PuOg, for the 1575, 3468 and 5679 w(th) sources, respectively. 


APPENDIX C 

236 Pu GAMMA PHOTON ACTIVITY 

The gamma photon activity due to 236 Pu contanimation in the Pu0 2 fuel 
was based on a 1.2 ppm impurity. 

The activities were calculated as follows; 

= 1.97 x IQ 13 x a "* >t x (abundance), 

sec-gm 23 ®Pu sec 


and 


dls 

sec-gm Pu02 


dls 

sec-gm 



x (am 236 Pu) 
(gm Pu0 2 ) 


The results of the calculations are tabulated as: 


(9) (9) 

Gamma Energy Abundance' 

(keV) , (% of isotope decay) 


48 

3.1 X 10~ 2 

no 

1.2 x 10‘ 2 

165 

6.6 x 10‘ 4 

520 

1.7 x 10" 4 

570 

1.0 x 10" 4 


. -4 

645 

2.4 x 10 


photons photons 

sec-gm ZiS<J Pu sec-gm Pu0 2 


6. 

107 

X 

10 9 

6, 

,46 

X 

10 3 

2. 

364 

X 

10 9 

2, 

,50 

X 

10 3 

13 

.0 

X 

10 7 

13, 

,7b 

X 

IQ 1 

3. 

35 

X 

10 7 

3, 

.54 

X 

10 1 

1. 

97 

X 

10 7 

2 

.08 

X 

10 1 

4. 

728 

X 

10 7 

5 

.00 

X 

10 1 
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The expression, e , was calculated to obtain the decay factor as: 


Time 0 1 5 10 18 

(years) 

Decay 1,0 0.784 0,296 0,087 0,013 

Factor 


APPENDIX D 

236 pu DAUGHTER NUCLIDE GAMMA PHOTON ACTIVITY 

The activity of the 236 Pu daughter nuclides was determined from the growth 
of 228 Th, using the Bateman equation. The decay chain for this family is: 


60m 0(66%) 3m 0 


0 . 3 M s 


The gamma photon activity is controlled by the growth of the 1.91 year 
228 Th. It may be calculated from the Bateman expression where 


236 Pu — 232u _JL_* 228 Th . a 

™ TTiSy u T2y^ in 1.91y 


is represented as: 


The final nuclide activity may be written as 


\jN 0 = X c N Ao X A *B e ~ XAt + X A*B e + X A X B e ct 

_(X B -X A )(X C ~^) (X a -X b )(X c ^B) (X A -X c )(X b -X c ) 

x (yield) x (abundance) x (concentration) 


i 





j 


where 


N 53 number of atoms of daughter C , 

n A0 , 


m a 

M a = actual mass of 236 Pu in Pu02 , 

M A = gram atomic mass = 236 gm , 

N 0 = 0.6023 x 10 24 , Avogadro Number , 

X A = . 24316 year -1 , ? 36 Pu decay constant , 

Xq = 0.009365 year -1 , U decay constant , 

- 1 228 

X = 0.36474 year , Th decay constant, 

c 

22 8 

The activity of Th as a function time was calculated as: 


TIME (years) 


dls/sec-gm 2 ^6p u 


1 

5 

10 

18 


2.738 x I0 10 
3.2216 x 10 11 
5.6957 x 10 11 
7.1566 x 10 11 
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When applied to p..n # , 

•he activity of , he ^ ™y vary, 


228 

Th Nuclide Actlvity/gni-Pu 0 

2 


236pu p°ncentration 

(ppm) 


2,8815 x 1 o 4 

* *02] x IQ 4 

1-441 x 10 4 
0. 240 x 1 0 4 


Vears 


3-383 x 105 
2.256 x 10^ 
1 -692 x 10 s 
2-819 x 10 4 


3*981 x IQ 5 
3*987 x 10 s 
2.990 x 10^ 
4 .984 x 1 q4 


U\- : \ ■ 


i 


I 

I 

I 

I 

I 

I 

i 
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APPENDIX E 

18 0(a ,n) 21 Ne GAMMA PHOTON ACTIVITY 

Tlw * 8 O(0< ,n) 2 *Ne reaction yields gamma photons as well as neutrons. If 
the total neutron emission rate for 238 Pu is assumed as 1 ,9 x 10^ n/sec-gm 
ZaaputH)^ 1,24 x 1Q 4 ,_Q.26 x IQ 4 and 0.4 x 10 4 n/seo-gm 238 Pu are 

the emission rate components for the * B Q(a,n) 2 *Ne, spontaneous iission 
and (a,n) impurity interactions, respectively.^ 37 ^ The photon activity for 
18 0 reactions may b© obtained as 

photons B 1.24 x 10^ n x photons /neutron emitted, 

sec-gm 238p u sec-gm 238p u 


and thus 


photons a. photon x 0.714 

sec-gm Pu0 2 sec-gm 23B Pu 


Gamma photon activities were calculated as 


Gamma Photon 
Energy (MeV) 


Photons emitted/ 
neutron emission M 


Photons 


Photons 
im PuO?. 





i 

1 

I 

These photon emission rates may be corrected for the decrease in 238 Pu 
I activity, using the decay factor e" *007734 year^t w hich yields 

I 

TIME 

| (years) Q 1 5 10 la 

I Decay 

| Factor 1.0 0.992 0.962 0.926 0.870 


I 

! 

1 

I 

I 

i 

I 
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APPENDIX F 

EXAMPLE FLUX CALCULATIONS 


An example gamma photon flux and an example neutron flux calculation 
are presented in this appendix. Both calculations are for a point located 
at a radial distance of 100 cm. from a 1575 w(th) capsule. Normal ( 18 0 
undepleted) oxygen is assumed in each case, 

1. The gamma photon flux for a 5 year old Pu0 2 fuel containing a 236 Pu 
impurity of 1.6 ppm is obtained as follows: 


From Table XVIII, the 'Decay of 
Isotopes and Daughters 'flux 

From Table XK the ' 18 0(a,N) 21 Ne 
flux 

From Table XX the ,236 Pu and 
Daughters' Flux x (1.6/1. 2) 


2,000 y/cm -sec 
2 

44 y/cm -sec 
5,181 y/cm 2 -sec 


Total flux 


8,952 y/cm -sec 


An estimate of the dose corresponding to this flux may be obtained, 
using Table 4-1, as 


* 17.54 x 0.175 + 17.54 x 0.825 (1.6/1. 2) 

= 22.4 mre m/hour 

2 Q c 

The above calculation assumes that only Pu change with age is significant. 
The oxygen flux yield may be corrected (x 0.962) using the table given in 
Appendix E. This calculation may be corrected for other capsule power 
values by interpolation, as in Figure 2-X. 
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2, The neutron flux for PuO^ fuel containing a 100 ppm boron (only) 
impurity is obtained as follows: 

From Table 3 -VII the 

18 0(«,n) 21 Ne + Fission Neutron' ® 427,7 n/cm 2 -sec 

flux 

From Table 3 -VII the 

'Fission Neutron' flux - 78, 5 n/cm. 2 -sec 

Thus,' 1 0(a,n) 21 Ne' flux » 349,2 n/cm^-sec 

Appendix Table K-ILgives the boron yield (maximum) as 46,9 n/sec-gm 
Pu per ppm or 4690 n/sec-gm Pu per 100 ppm, From Equation (3.5) the 
boron impurity flux is determined as 

= 349.2 x 4,690/9,882 n/cm 2 -sec 

“ 156.7 n/cm 2 -sec 

The total flux is simply equal to 

165,7 + 427 = 593.4 n/cm 2 -sec 

An estimate of the dose rate corresponding to this flux may be obtained 
by scale-factoring the 52 , 4 mrem/hour dose rate given in the text of 
Section 4 as 

- 52.4 x (593.4/427.7) = 72.7 mrem/hour. 

The above may be obtained for. other power values by repeating the cal- 
culations for 3468 and 5679 w(th) and interpolating at the desired power. 



APPENDIX G 

THE DYNAMICS OF THE 18 0(a,n) 21 Ne REACTION 

The (<*,n) reaction may be considered as an inelastic collision between two 
initial systems and the reaction written as 

IQ 21 

0*n+ Ne + Q 
8 10 

The total. kinetic energy of the interacting system corresponds to the kinetic 
energy of the alpha particle 

E a = E f - Q 


where 

Q is the energy equivalent of the difference between the masses of the 
initial reactants and product involved. 

It can be shown that the total kinetic energy before the collision in the cen- 
ter of mass coordinate system' is the sum of the. kinetic energy of the incident 
alpha particle and oxygen target nucleus. The reaction masses and energies 
may be defined as follows: 

M a = the rest mass of the alpha particle, 

M 0 = the rest mass of the oxygen atom, 

Mn = the rest mass of the neutron, 

Mn 6 “ the rest mass of the neon atom, 

Ea = the kinetic energy of the alpha particle in the laboratory system , 


I 


I 

I 

1 


l , 
> 

i ! ; 


f :■ 


0 ' 
p 




<r; 

A * 


'7 


i 


t 

l 


E a > « the kinetic energy of the alpha particle in the center of mass 
coordinate system, 

E n = the kinetic energy of the neutron in the laboratory system of 
coordinates , 

E n ' » the kinetic energy of the neutron in the center of mass coordin- 
ate system. 


The total kinetic energy prior to collision may be written as 


E 


i 


t 1 

p M 0 

t, > 

“ \M 0 + Ma 


the total energy of the incident alpha particle relative to the target nucleus . 


The total energy of the incident alpha particle in the center of mass system 
is 


Ed 


[ Mq 

\M 0 + M a 


The energy of the emitted neutron in the center of mass coordinate system 
is obtained as 



M NE ^ 


Q + 

fl- M a 1 

E« 

E n “ | 

i M NE + M NJ 


[ m n + m ne j 


which reduces to 


E 


t 

n 


m NE I 

l m n + m ne I 


IQ + M p Eg 

Mq + 


m NE 
+ M 



(Q + E a ') 


The emitted neutron onorgy in the laboratory system of coordinates is E n 


where 


. m N M a. s d + HjjMg 

fi n"*n r (Mq+Mo ,) 2 (M 0 + M ft ) 


2E<* / 2E n ' cos 0cm 


This expression relates the neutron energy in the laboratory system,. E n , 
to the incident alpha particle energy, Eg and the angle of scatter, 0cm, 
in the center of mass coordinate system. 

The kinetic energy of the emitted neutron is determined from the Q value 
of the reaction, the incident alpha particle energy, Eg, and the angle 
of scatter in the center of mass coordinate system as 

E n = a (Eg) + b(E a ) cos 0cm 


where 


+ M 


Q + M 


+ M 


+ Mjt M(* 


(M 0 + Mot) 


(M 0 +M a ) \U a 


M N +M NE\ M o +M a 


If the probability of emission into an angle, 0cm, in the center of mass co- 
ordinate system, is isotropic, then p(0cm) (sin 0cm)/2, is the probability 
of scatter such that 


[•Yr 





i 



p(0cm) dflbm » 



o 

a In 8cm 


2 


dQcm ® 



1/2 d(eos 0cm) 



l/2dM=l 


where 


djU” d(cos 0cm) . 


Rewritting 


sin 0cm « (1- cos^Scm) 1 / 2 

M 


= /l-/k n - a\ 2 ] 1/2 = /b 2 -E n 2 + 2E n a - a 2 ^ 1/2 


the probability of scatter can be rewritten 



d E n =l 


This equation gives the probability of scatter in the center of mass coordin- 
ate system in terms of a(Ea), b(E a ), and E n which are parameters of the 
alpha particle and emitted neutron energy in the laboratory system of coordin- 
ates. It integrates to unity for particles scattered over all angles of emission. 
The equation expresses a uniform distribution of scattered neutrons in the 
laboratory coordinate system. 
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APPENDIX H 

ALPHA PARTICLE ENERGY LOSS PER UNIT PATH LENGTH 

238 

Alpha particles are emitted during the decay of * 00 Pu with an energy of 
5.50 MeV. They eventually lose this energy by collision with the oxygen 
and plutonium atoms present in the Pu0 2 fuel. The loss in energy, per 
unit of path length traveled by the alpha particle, may be determined from 
the expression: 



dE a 


\ 




3.8117 x 10 3 


J5L Niln [5.4823 x 10 2 E^ 

E 0 [ tl J 


where 

1^ = the mean excitation energy of the material species, i. 

The mean excitation energy was assumed to be 91eV and 893eV for the 
oxygen and plutonium isotopes . 

E<* = the energy of the alpha particle (MeV) 

Zi = the atomic number of the scattering nucleus i, 

3 

= the atomic density of the scattering nucleus (atoms/cm 0 ). 

A more detailed discussion of this subject is given in the excellent article 
by Fano ( 36 ). 


APPENDIX I 

18 0(a, n) 21 Ne NEUTRON YIELD 

Tha 5,5 MeV alpha particles emitted by the radioactive decay of plutonium 
238 interact with AO Q by a inelastic collision to produce neutrons. The 
alpha partial© in this (a,n) reaction will have a maximum energy of 5.50 
MeV and a possible minimum threshold energy 

Eth “ "(1 + mg ) Q 
m u 


Q is the energy equivalent of the difference in masses of reactants and 
products of the reaction. This term is discussed in more detail in Appendix 

G. 


The incident alpha particle with energy E th * E th * 5.5 MeV collides with 
18 0 and forms the compound nucleus 22 Ne. The 22 Ne emits a neutron and 
decays to the excited and ground state levels of 2 *Ne. The yield and Q 
values from this reaction were given in Section 3.3 as: 55% to the ground 
state (Q = -0.70 MeV), 35% to the first excited state (Q = -1.05 MeV), 
10% to the second excited states (Q = -2.43 MeV), and 1% to the third 
excited state (Q = -3.54 MeV). 


The number of neutrons produced in each energy group j and the correspond- 
ing Q value of the reaction was calculated according to the relationship: 


Nj(Q) 


r r* (Ej) 

“J a ( E a^ E a I 

J Ea -Wj+1) 


1/2 du dE, 


a 


3.8117 xIO 3 |" 16log e 5 ' .4823 xlO 2 Eaj+ 941og e 5. 4823x1 0 2 Eg j 


r 
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where 

CT (*•£*) " the cross section for the formation of the 22 Nq ( 39 \ 

E a = the energy of the alpha particle in MeV. 

M(Bj+i)/ M(Ej) are the lower and uppor values of cos0 cm from which 
the energy group j can be produced . 

The number of neutrons produced in each energy group J , by the decay of 

the compound nucleus through 22 Ne all possible decay channels, was 
calculated as 

Nj = 0.55 Nj (-0.70) + 0.35 Nj (-1.05)+ 0.10 Nj (-2.43) + 0.01 Nj (-3.54) 
where 

Nj (-0.70), Nj (-1 .05) , Nj (-2 .43) , and Nj (-3. 54) are the number of 
neutrons produced from the particular Q value reaction. 

The sum, SNj is determined as the total number of neutrons produced in the 
0(a,n) Ne reaction. The fraction Nj/pNj as calculated for an arbitrary 
10,000 groups j, is plotted in Figure I-I. 

The equation relating the neutron energy, E n , the incident alpha particle 
energy, E a in the laboratory system of coordinates, to the angle of scatter, 

0 C m in the center of mass coordinate system is 


E n = a(Eflr) + b(E a ) cos Ocm 


This equation has been derived and discussed in Appendix G. 

To illustrate the relationship between E n , E a and 0 cm (angle of scatter), 
Figure I-II was tabulated for the following conditions: 

Q » -0.70 MeV; the Q value for the inelastic collision which will 
yield the ground state of ^Ne. 

M NE (mass of isotope) « 21 .00000 a ,m .u. , 

Mjj (mass of incident neutron) - 1 .008665 a .m .u . , 

Mo (mass of target nucleus ^®0) = 17 .999160 a .m .u. , 

Ma (mass of alpha particle) = 4 .001506 a .m.u. where a ,m.u. is the 
standard atomic mass unit. 
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NEUTRON FRACTIONAL YIELD AS A FUNCTION OF 






APPENDIX J 


INTEGRATION OF 18 0 (a, n) 21 Ne SPECTRAL YIELD 


where 


and 



t * 

Nj ® 0.55 Nj <Q = -0.70 MeV) + 0.35 Nj (Q = -1.05 MeV) 

+ 0.10 Nj (Q * -2.43 MeV) + 0.01 Nj (Q « -3.54 MeV) (3.6) 


£ 


18 21 

Nj = the total number of neutrons produced by the 0(<*,n) Ne 
reaction. 


N(Ej) « the numberxif neutrons in the laboratory system, between 
El and Ejj, the lower and upper energy limit of the interval. 


The number of neutrons in each group j, for each Q value, was numerically 
integrated according to the relationship: 



where 

M(E) = cosine of the scattering angle in the center of mass 
coordinate system, 


;i07 


■ the number of neutrons in each group j in the laboratory 
system of coordinates (j « 10,000 was used in the cal- 
culations) , 

1 8 

0(E/J « the compound (<*, n) nuclear cross-section for the O 

(39) 

raactien corrected to the laboratory system of co- 
ordinates, 

E a ■ the range of the alpha particle energy appropriate to group 
j and each value of Q. 

Since isotropic scattering is assumed in the center of mass coordinate 
system, the integral in equation (3.7) 


Hi Ej) n (Ej) 

I 1/2 d u = j 1/2 d(cos0Qjj^) ® j 1/2 sin^j^ ^®cin * 

•i(E J+1 ) iOj+i) J.(b, +1 ) 

corresponds to the neutron distribution, in the C. M. coordinate system, 
contributing to group j in the laboratory coordinate system . 

tithe compound nuclear cross-sections cr(Eg) are plotted as a function of 

(39) 

incident alpha particle energy, Shapiro suggests that although the dis- 
tribution will have the correct shape, the numerical values of the compound 
nuclear cross sections may be low. To correct for this fact the fraction 
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4 238 

was calculated and multiplied by 1.22 x Ur neutrons/sec-gm Pu, the 
experimentally measured yield of the 18 0(a ,n) 21 Ne inelastic neutron re- 
action rate in PuOg (the normalizing value of the total yield from the 
1 8 Q(a,n) 2 ^-Ne reaction). 

The nuclear cross sections for the. formation of the compound nucleus 2 qNo 
were determined from the tabulated quantities of Table J-I , with g = 4 and-.. 
v Q / B = 1 -0 as given by Shapiro^ 39 ) . The cross sections calculated by this 
method are normally given in terms of the alDha particle energies in the 
center of mass coordinate system. Table J-I . gives the alpha particle 
energies, converted fi. m the C. M. to the laboratory system and the corres- 
ponding calculated nuclear cross sections . 

Although the data of Baird and Willard^ 33 ) is an excellent source of nuclear 

1 ft 

cross sections for the 0(a,n) reaction, the lowest alpha particle energy 
is 2.25 MeV and thus the data is difficult to interpret for averaged cross, 
sections near the resonance psaks; data for alpha particle energies down 
to at-least-1.0 MeV is desirable-. 

If the relationship between neutron energy .group, j , alpha particle energy, E a 
and the angle of scatter, cos6 n m « p , in Appen dix G , is written as 

M(Ej) -OiEj + 0 , 

then the following combination of integration limits apply for each neutron 
energy group j and vab <4 Q, where, m(Bj) = (+1)* (-1) anc * (** “ 1 and <+ *) 
are defined as cases A, B and C, respectively: 

M(Bj) - M(Ej+i) 

A-A ■ (1 - l)/2; no neutrons in group j 

A-B » (1 -(- l)/2 =1; all neutrons are in group J 
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A-C » (1 -(«Ej +1 + /3) )/2; some of the neutrons are in group J 

B-A cannot occur since, #*(Ej) >/fi(£j + i) 

B-B « (-1 -(-1 .01))/2; no neutrons in group J 

B-C same as B-A 

C-A same as B-A 

C-B = (oe Ej + /3 + l)/2; some neutrons are ingroup j 

C-C « <*(Ej -Ej+i)/2; some neutrons are in group J 
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TABLE J-I 

» 

NUCLEAR CROSS SECTIONS FOR FORMATION 
OF 18 0 (<*, n) 21 Na COMPOUND NUCLEUS* 39 * 


Eg (MeV) g(Ea) (10~ 24 bams) 


L.222 

2.52 x L0 -3 

1.834 

10.69 x 10~ 3 

2.445 

17.32 x. 10" 3 

3.056 

48.36 x 10 -3 

3.667 

23.7-5 x 10“ 3 

4.278 

166.04 xJUf 3 

4.890 

189.30 x 10” 3 

5.501 

198.13 x 10‘ 3 
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APPENDIX K 

CONTAMINANT LIGHT ELEMENT (ft,n) NEUTRON YIELD 

Low atomic number elements present as contaminants M, which ara. 

uniformly mixed with the plutonium isotopes in the commerical Pu 02 

fuel , produce additonal neutrons by the (o< , n) reaction with the impurity 

. . . ( 3 , 53 , 54 ) 

given element . 

The neutron yield Ny, for the X(« , n)Y reaction with the given Impurity 
element X, where X is uniformly mixed with the Pu02 , may be obtained 
as 


Ny = KXNpyJEjj 


Sx N x 

SpuNpu + S X N X 


where 


the yield is the number of neutrons produced per gm of Pu per million 
parts by weight of the contaminant element, X, 

K = the 238p u fraction of plutonium, 

A = the ^^®Pu disintegration constant, 

N x = the number of atoms of element X, 

Np U - t:..o- number of atoms of the plutonium present, 

T x = the thick target yield for the element X (yield produced by the 
alpha particle incident on the thick target composed of the 
element X) . 

S x = the atomic stopping power for the alpha particle incident on 
material species, x. 


Neutron yields were calculated using the thick target yields and atomic 
stopping powers data given in Table K-I. 
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The (tt,n) neutron yield Ny was calculated using the approximate rela 
tionship 


N y - KXT x 


1 

\l*s] 

l A *J 



» 1.4807 X 10 8 K 



where 


N 0 « 6.024 x 10^3, Avogadro number. 

A x = the atomic mass of the element x. 



If minimum and maximum values of the ratio of the atomic stopping powers 
are multiplied by the minimum. and maximum thick target yield, respectively, 
Table K— II yields result. 
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TABLE K-I 

DATA FOR LIGHT ELEMENT (a ,n) YIELD 


Element 

Ratio of Atomic Stopping 

Thick Target Yields (in units 


Powers (S x /Spu) 

of 10" 6 

neutrons /alpha) 


Sx/Spu^* 

8x/Sp U (52) 

^(54) 

Ix (53) 

Li 

.1038 

.13 

2.4 

2>4 (26)»* 

Be 

.1320 

.17 

72. 

84.4 ± 0.9 

B 

.1579 

.20 

21. 

19.6 * 0.2 

C 

.1821 

.23 

.09 

0.113 ± .015 

N 

.2047 

- 

.01 — 

mm 

0 

.2260 

.29 

.061 

.068 ± .011 

F 

.2462 

.31 

10.4 

11.6 ± .2 

Na 

.2837 

.34 

1.3 

1 , 5 (28)** 

Mg 

.3013 

.36 

1.2 

1.33 ± .04 

A1 

.3180 

.37 

.64 

.76 ± .03 

Si 

.3342 

• 42 

.15 

.168 ± .020 




where Z = atomic number . 


** these uncertainties were given as +25%, -40%. 
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TABLE K— II 

LIGHT ELEMENT (a , n) YIELD 
(rveutrous/see-gm Pu, per ppm impurity) 


Element 

Maximum Yield* 

Minimum Yield* 

Yield 

Li 

5.5 

4.3 

4.6 

Be 

191.0 

126.5 

133. 

B 

44.9 

35.1 

41. 

C 

.3 

.2 

.2 

N 

- 

mm 

- 

0 

.2 

A 

.1 

P 

22.7 

15.2 

18.0 

Na 

2.7 - 

1 A 

2.2 

Mg 

2.4 

1.8 

2.1 

A1 

1.3 

.9 

1.0 

Si 

.3 

.2 

.2 


* The maximum and minimum yiei- 5 ' . e not observed yields, but calculated 
using all possible thick target yields and stopping powers given in the 
literature . The values will include uncertainties in the observed thick 
target yields. whenuaecLby the reader. .. ..... 
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